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' Abstract 

I This is a contribution to the Proceedings of the 7th Mathematical Physics Meeting: Summer 

^S) . School and Conference on Modern Mathematical Physics, held in Belgrade 09 - 19 September 2012. 

> ^ ' We give an easily accessible introduction to the ADM decomposition of the curvature components. 

, After that we review the basic problems associated with attempts of constructing a viable massive 

gravity theory. And finally, we present the metric formulations of ghost-free massive gravity models, 
and comment on existence problem of the matrix square root. 



(N 



O 



1 Introduction 



General relativity is a very successful theory. We are not aware of any experimental contradictions to its 
predictions whenever we are able to test it directly. However, at the largest scales, it fails unless we are 
willing to introduce huge amounts of Dark Matter and Dark Energy. The Dark Matter must be of some 
yet unknown mysterious nature. And the Dark Energy might be a mere cosmological constant, but of 
an overwhelmingly unnatural value. At the end of the day, it pushes us towards reconsidering the basic 
J>. ', foundations of the theory of gravitational interactions. 

' Modifying the General Relativity, even if just by giving the graviton a mass, is not an easy game to 

I play T. And it is only very recently that we have got a decent hope [2, 3 of obtaining a viable model 

of massive gravity. In this contribution we review some non-perturbative aspects of this theory from the 
viewpoint of the Hamiltonian ADM analysis. 

(N 
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^ ! 2 Review of the ADM formalism 

In the ADM approach, we use the following (3 + l)-decomposition of the metric [4] 
X ■ ds^ = - {N^ - NiN') dt^ + 2Nidx'dt + j.jdx'dx^ (1) 



where N and Ni are commonly called lapse and shifts respectively, and the spatial indices are raised 
and lowered by the spatial metric 7. Equivalently, goo = ^ {N^ ~ NiN'-), goi — Ni and gij — -fij. And 
it is fairly straightforward to check that the inverse metric components are g'^*' — g^^ = and 

gij — ^ij _ iLj^. Comparing this result with the Cramer's rule for 5'^" we get the useful relation for 
determinants, y/—g = Ny/j. 

Let us also introduce a unit normal vector n'^ = (^jf,^^^ which corresponds to the one- form 

= (^—N, ^ nullifying the tangent vectors. Generically, these objects are not covariantly constant, 
and one can define the extrinsic curvatures of the spatial slices t — const: 

if,, ^ - V. n, = r^n^ = -NT^ (2) 

with covariant derivatives ^ ^riy = S^n^ — V^^ria and ChristofFel symbols 



of the metric g. 
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(3) (3) 

We easily get Tijk — T ijk where r denotes thcl Christoffel symbols for the metric 7, and Foy 
i {diNj + djNi — jij) which gives 



1 /(3) _ (3) 
27V 



K,^ = -iVrO^. = — ( ^,N, + ^,N, - 7,, ) (3) 



for the extrinsic curvatures ^ where the three-dimensional covariant derivatives are taken with respect 
to 7. Incidentally, we note that 



— n X . (4) 
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In order to analyse a gravitational theory, we need to compute the curvature components 



R^.cp = 5.r^, - 9^r^j, + rjj.r^, - r^/s, (5) 

in terms of the ADM metric decomposition ([l}. For the Levi-Civita connection we have the symmetry 
properties Rap^u = —R/Safj.i' = —Rapu^i = Rnvap, and therefore need to know only three types of compo- 
nents Riju , Roijk and Roioj . The calculations are rarely presented in an explicit and detailed form. One 
can either use the geometric meaning of the curvature tensor and play with parallelly transporting normal 
and tangent vectors, or one can follow more direct derivations in a somewhat involved algebraic disguise 
in mathematical textbooks. We want to present a straightforward, almost brute force, computation. Our 
main point is that it can also become rather simple if cleverly done. 

2.1 Connection components 

It is reasonable to calculate the connection components first. Anyway, they might be needed for the 
matter part of the action. We will thoroughly trade 7 for the extrinsic curvatures using ([3]). In the 
following we need: 
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N^K,,) , (9) 
-N[l''-^;^jKk, + V,N\ (10) 

(11) 

Kjk- (12) 



2.2 Riemann tensor components 

Using the formulae (|6]) - (fT2)) . we readily obtain the curvature tensor ([5]). First, 

Rijki = gipdk^ij - gipdiT^j + r^fepFf^. - TtipT^j 



I 1 \ /(3) AT™ \ 1 / (3) 

-N,dk ( -K,ij +j,„A ir]l + —KjA - -K,i ( ~NK,k + SJkN, 



(3) /(3) ATI 

T,km[ TT,+—Kvj]-{k^l) 



(3) 

= Rijki + KikKji — KuKjk (13) 



Then it would be easier to calculate 



1 A''*^ 
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instead of Roiaj. Indeed, we get 

(3) 



(3) (3) 

\J jKki- \JkKji (15) 



Finally, after just a little bit of very simple algebra, we find 



(3) (3) (3) (3) 

Using (jl5l) and (1141) we trivially transform it to a more symmetric form 

1 / • (3) (3) \ 

n'^n'^i?^,,, = — Ay- + \j, \j,N + NKi^Kkj - Cie^K,, 1 (16) 

where Ciej^Kij = N'^dkKij +KikdjN'' + KjkdiN'', and the partial derivatives can be substituted by the 
covariant ones. 

2.3 Scalar curvature and Einstein-Hilbert action 

The formulae - ([T5)) represent the full Riemann tensor for the Levi-Civita connection. Now the 
Ricci tensor Rf^^, = R°' ^ai/ components n^Rf^i and n^^'nyRf^i, can be easily obtained. We skip this 
point and find the Ricci scalar R = g'^'^R^i, directly. Using and symmetry properties of the Riemann 
tensor, and the obvious relation 'f^^ Kij — OqKI + K^^'jij together with ([3]) we get 

R = g^-'g^'^R^^.p = f''l''R^,kl - 2n^n''f^R^,,, 

(3) 2 ■ ■ • 4 --(3) iV-''(3) 2 (3) 

(3) . ■ 2 • • A^J' (3) . 2 (3) 

= R+ K^^K,, + KIK^ -Kl + 2— ^ ,Kl - ^ A A^ (17) 

As a final step, let us put the Einstein-Hilbert density ^—gR — ^/jNR into a convenient form. We 
use the relation ~ -^Y^iij and exclude 7 by virtue of ([3]): 

/(3) . A . • . (3) 

V^R = f i? + K'^K,, - KlXn - 2^g - 2^^N 

(3) 

where -27=^ Vm (^X) -29o (Vt^^O (i^l^^)- 

Neglecting total time derivative and covariant divergence terms, we obtain the Einsten-Hilbert action 
in the ADM formalism: 



S = j dtdPx^N {^R + K'^K,, - KIK] 



(18) 



The action (|18p would be enough for us as the kinetic part of the massive gravity models. However, as we 
have computed all the curvature components, these results can be used for more complicated modified 
gravity models, too. Note also that the curvature tensor in theories with non-metric connection can be 
viewed as an (exact at second order) variation of the Riemannian one with respect to the connection 
being varied from its Levi-Civita value by non-metricity and contortion tensors. 



3 Massive gravity and its problems 

Now we come to massive deformations of general relativity. Note that only the six 7^ variables are 
dynamical in the action (fT8|) . The lapse and shifts play the role of Lagrange multipliers which enforce 
four first-class constraints corresponding to the diffeomorphism invariance of the theory. The total number 
of physical degrees of freedom is two. A generic mass term would break the gauge invariance and produce 
a model with six independent degrees of freedom. However, a massive spin-two field should posess only 
five of them. Therefore, we have an extra scalar. 
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That this scalar might be problematic, can be seen already at the level of quadratic perturbations 
around Minkowski space. For linearised theory, the Einstein-Hilbert action acquires the form of 

-l{do.h^.){d''hn + lid'-h.^Wh^J - l{d^h'^^){d,h^p) + ^id.hZKd^hf'p) 

which can be studied in the standard parametrisation of /loo — '^4'^ hoi — dib + Si with diSi = 0, and 
h,j = 2i'd,j + 2(92 + d,Vj + djv, + h[J^^ with d,v, 0, d.hlj^^ = and h\f^^ = 0: 

- J(5a^^r^)(5"/i!r') + I (dj {v^ - s,)f - 6V? + 2(9,7A)2 + 4V^ A - 6 + a) 

where we have neglected all total derivative terms. One can see that there are two gauge invariant variables 
h\j in helicity-two (transverse traceless) sector, two gauge invariant variables Vi — Si in helicity-one 
sector, and two gauge invariant variables, '0 and (j^ — b + a, in helicity-zero sector. Of course, these are just 
the standard gauge invariant variables of the cosmological perturbation theory in the limit of vanishing 
Hubble constant. Only the helicity-two modes are physical, the others are constrained. And the latter is 
extremely good because otherwise we would have got severe problems with the wrong-sign kinetric term 
of ip. Unfortunately, we do indeed get them once the gauge invariance is broken by, say, a mass term. 

In general, we can think of two types of a mass term, h^^h^^ and h^h'^- Both of them contain /igo 
and /loi fields non-linearly (quadratically) which implies that the corresponding Einstein equations cease 
to provide constraints for the spatial sector. However, there is one particular combination of them, 

discovered by Fierz and Pauli [5], from which ft,§Q drops. The resulting constraint kills the unwanted sixth 
degree of freedom and gives a ghost-free massive gravity in the linear approximation. 

Unfortunately, even in the vanishing mass limit, the linear theory anyway contradicts observations [B] 
due to the scalar graviton which couples to dust modifying the effective gravitational constant, but not to 
radiation keeping the bending of light intact. It was later argued by Vainshtein [7] that non-linear effects 
will take over at small scales and restore the GR limit. But, almost at the same time, Boulware and 
Deser have shown [8] that the sixth degree of freedom comes back at the non-linear level reintroducing 
the ghost mode. And therefore, a stable theory of massive gravity is probably not possible at all. 



3.1 de Rham-Gabadadze-Tolley model 

Recently, a potentially ghost-free non-linear construction of massive gravity was found [2] in the decou- 
pling limit (m — >■ and Mpi — ^ oo with rn^Mpi fixed) by carefully getting rid off the sixth mode at every 
order of perturbation theory starting from the Fierz-Pauli term. The model can be resummed [3J H] to a 
deceptively simple potential 

V = In? (TrV5~V - s) (19) 

around any background metric /^j^. Up to the —3 term which is needed to avoid a linear in h = g — f 
contribution to the action, this is the first symmetric polynomial of eigenvalues of the square-root matrix 
\f(f^ ■ (Interestingly, it is noted in [TU] that this model was known to Wess and Zumino as early as 

in 1970.) The other two suitable potentials are the second ^Tr^/ <?^^/) — Tr g^^j^ and the third 
(^TrA/g-i/)^ - 3 (Tr 1/5-1/) Tr {\J g-^ + 2Tr {\Jg-^i^ symmetric polynomials. 



4 The Hamiltonian analysis 

In order to count the degrees of freedom non-perturbatively, one has to perform the Hamiltonian analysis 
of the action (fT8|) with the potential term —^NV given by (IT9l) in the case of the minimal dRGT model. 
We find the momenta tt*^ = — ^ (^Kj^j^^ — K^^), the primary constraints ttn = t:^. — 0, and the 
Hamiltonian 

/ ("^ ) \ 
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For a general potential, the lapse and shifts enter the Hamiltonian non-linearly, and therefore, the sec- 
ondary constraints (obtained, as usual, by commutation of the Hamiltonian with the primary constraints) 
directly determine their values instead of constraining the six spatial degrees of freedom. The major claim 
is that for the dRGT potentials there is a combination of the lapse and shifts which continues to serve 
as a Lagrange multiplier enforcing a constraint equation on the spatial sector. 

4.1 Hassan- Rosen proof 

Hassan and Rosen have found \T2\ a redefinition of shifts 

N' = {S)+ND}{j,n))n' 
such that, for the new shifts n', the square-root matrix takes the form 

1 f I n*\/0 



NVT^W [ -n'n^ j ^ I X^'{l,n) ) ^^^^ 

which means that the lapse assumes the role of the Lagrange multiplier and imposes a non-trivial con- 
straint on the spatial sector. One can check that the same is true for the higher ghost-free potentials 
because the first matrix in (j20p preserves its form under exponentiation to an integer positive power. 
A careful investigation shows [13] that the subsequent consistency conditions provide the second spatial 
sector constraint to form a non-degenerate pair of second class constraints on 7^^ , and finally, an equa- 
tion for determining the lapse. And the same is true |14j for bimetric versions (with an independent 
Einstein-Hilbert term for f^u)^ too. 



4.2 Auxiliary variables 

One can also approach the problem in a different manner. We can introduce a matrix $ of auxiliary 
fields with a constraint that "I>^ = N'^g^^f which amounts to 

NV = 2mH'; + < (ci>^$^ _ N'g'-^U,) 

and, varying with respect to $, integrate out the k-s: 

NV = + {<P-X iV2g-/„^) . 

In this method it is not necessary to calculate the matrix square root explicitly, and the procedure is 
standard and elementary. Unfortunately, the calculations are cumbersome. But one can easily find [15) 
a direction in the space of must-be unphysical variables {N, iV% $) along which there is no restriction 
at the level of secondary constraints which proves that the number of degrees of freedom is less than six. 
Note that, in this construction, the equations of motion automatically demand that a real square-root 
matrix \fg~^S does exist. 



5 Existence of the square root 

Non-pcrturbative existence of the real square-root is one of the very interesting issues about the dRGT 
mode [16] . In the cited paper the reader can find the necessary and sufficient condition for it: if there is 
a real negative eigenvalue of the j matrix, then it must be accompanied by even numbers of identical 
Jordan blocks. However, let us note that if f has a negative eigenvalue —A, then the corresponding 
eigenvector lies in the kernel of At; -I- / matrix. In other words, there is a degenerate linear combination of 
g and / with positive coefficients. Even though it is anyway undesirable to couple matter to a non-trivial 
linear combination of metrics, one might consider imposing a condition of non-degeneracy of all positive 
linear combinations. On the other hand, in full quantum gravity it might appear to be not so much 
non-sensical to consider complex potentials. 

We would point out that there is also a uniqueness problem for the squrc root. Even the 2x2 unit 

matrix has an infinity of real square roots, for example: ^ that 

this one has zero trace which shifts the vacuum energy with respect to the trivial square root. It might 
play an important role in a fully quantum regime. 



5 



In conclusion, let us mention that there are lots of other problems to take care of. Some of remaining 
five modes may become ghosts in certain regimes, and there can be strong coupling issues in some phys- 
ically relevant solutions. And there are interesting problems concerning multimetric theories. Moreover, 
we have to accept some form of acausality in our physical world since the superliminal propagation was 
shown to be quite generic in the models of dRGT type ^10) . However, it is very remarkable that, in the 
first place, we do have a reasonable candidate for a viable massive, or bimetric, gravity model. 
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